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Abstract 

We study ab initio computations of the interaction of Lithium with a strong laser field. Numerical 
solutions of the time-dependent fully-correlated three-particle Schrodinger equation restricted to 
the one-dimensional soft-core approximation are presented. Our results show a clear transition 
from non-sequential to sequential double ionization for increasing intensities. Non sequential double 
ionization is found to be sensitive to the spin configuration of the ionized pair. This asymmetry, 
also found in experiments of photoionization of Li with synchrotron radiation, shows the evidence 
of the influence of the exclusion principle in the underlying rescattering mechanism. 
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Photoionization of atoms by short-pulsed intense laser radiation constitutes an extraor- 
dinary playground to test quantum mechanics beyond the perturbative limit. The earliest 
experiments on this subject already showed many new features unexplained by the standard 
approaches of that time, and renewed the interest on alternative non-perturbative theoretical 
developments. Among them, the ab initio numerical integration of the Schrodinger equation 
provided a fundamental tool for the profound understanding of the dynamics of laser-matter 
interaction. Always limited by the current state-of-the-art of the computing hardware, the 
numerical codes first targeted the problem of single electron ionization in one dimension 
using a soft-core potentia! Q Q tffl Q. The Motion to ID has heen considered for 
more than ten years a reasonable strategy to get a first insight into the dynamics of intense 
field ionization. In particular, many aspects of the strong field phenomena can be described 
qualitatively at this level (among them, above threshold ionization, and the basic general 
features of the harmonic spectrum). More recently, one-electron three-dimensional codes 
have helped to unravel more subtile phenomena as non-dipolar effects [7J. However, the 
complexity of the numerical task grows exponentially with the number of particles, as new 
dimensions have to be added. This introduces a rather serious technical limit to the compu- 
tation of the full 3D dynamics of more than one particle. Nowadays, the exact integration 
of the 3D Schrodinger equation can be accomplished only for the case of He interacting with 
linearly polarized electromagnetic fields, employing a extraordinary amount of computing 
resources 0, 0). In these circumstances, the dimensional reduction of the many-particle 
problem continues as a fundamental tool. For instance, the ID approach to the laser-He 
interaction j^, |j| is still employed as the most common technique to tackle the two-particle 
problem. For three electron problem, the dimensional reduction appears almost mandatory. 
Note that other traditional approaches (density functional theories) have not straightfor- 
ward applications in the limit of small number of particles Q| , which are highly correlated. 
Quantum correlations, therefore, play a fundamental role in the dynamics of few particles. 
The advantage of the dimensional reduction is to allow ab initio numerical calculations that 
include completely these correlations. 

Accordingly to the underlying mechanism, the double photoionization of Helium can be 
cataloged as sequential or non sequential. In the first case, both electrons ionized indepen- 
dently by photon absorption from the electromagnetic field. In contrast, the non-sequential 
ionization reveals a more subtile dynamics, in which the second electron is ionized via scat- 
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tering with the first ^} 12, Q]. One signature of the relevance of quantum correlation in 
this later process consists in the sensitivity to the particular form of entanglement of the 



ionizing pair. As reported in 
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la ], the rescattering process is less effective when the 



two-electron wavefunction is antisymmetrized in the orbital part (orthoHelium) rather than 
in the spin part (paraHelium) . In the three particle problem, i.e. Lithium, correlations 
appear more intrincated involving non-separable orbital and spin antisymmetries. Exper- 
imental works on the double and triple ionization in Lithium has been published recently 
for synchrotron radiation and ion or electron collision 3, 18, Q,!^. On the other hand, 
previous theoretical treatments include the high photon energy limit [22], approximated 
half-collision models 24 , 2^ and, very recently, close-coupling grid calculations in the weak 
field limit 2(J- 

To our knowledge, the problem of Li photoionization in strong laser fields using ab initio 
numerical calculations of the three-particle problem has not been previously addressed. Of 
course, the full 3D problem falls well beyond present and near future computing capabilities. 
However, the problem in reduced dimensionality (ID for each particle) can be addressed with 
a medium-size computer. We, therefore, present in this paper the first results of this type 
of calculations that consider fully correlated electrons. 

In the limit of very high photon energies, "shake off" has been determined to be the 
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However, in the case 



main mechanism for double and triple ionization of Lithium 
of photon energies below some hundreds of eV, a different mechanism has been proposed 
(20, S, Q] • In this case the electromagnetic field ionizes one or two electrons from the inner 
K-shell which, in the way out, ionize one of the remaining electrons. This viewpoint seems to 
be confirmed experimentally by Wehlitz et al [17], with synchrotron radiation. In addition, 



the same experiment suggests through the comparison with photoionization of He, that the 
double ionization of Li is not equally efficient for the different spin configurations of the 
ionized pair of electrons. The present study confirms this aspect and gives a fundamental 
description in terms of the inhibition of — > 2e~ scattering due to the Pauli's exclusion 



principle. As stated previously 



3,0, 



the symmetric character of the spatial wavefunction 



with respect to exchange of particles, can inhibit non-sequential double ionization. 

We construct the three-particle hamiltonian in reduced dimensionality by extension of 
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the previous models for Hydrogen and Helium atoms (in a.u.): 

3 \ ^ 1 (1) 

where a and b are the parameters of the soft potential. This form of hamiltonian commutes 
with the symmetry operators and, therefore, the symmetry of the wavefunction remains as 
a constant of motion. Initially, we will assume the atom in its ground state 2 Si/ 2 , therefore 
the wavefunction at any time has this symmetry and may be expressed as 

(z 1 ,z 2 ,z 3 ,t) oc a(l)a(2)/3(3)0i 2 (zi,Z2,Z3,*) 
+ a(l)/3(2)a(3)4> 13 (z 1 ,z 2 ,z 3 ,t) 
+ P(l)a(2)a(3)(f> 23 (z 1 ,z 2 ,z 3 ,t) (2) 

The spin part is the combination of three single electron spin functions, in our case a(i) = 
| I) and (3(i) = | |). The orbital functions 4>ij(zi, z 2 , z 3 ,t) are antisymmetric under the 
permutation i j. Note that we have written (J2J) in such a way that the different terms 
in the summation have orthogonal spin states. This form is particularly useful with the 
non-relativistic hamiltonian ([TJ since the spin state is a constant of motion and, therefore, 
every term in the summation evolves independently from the others. Moreover, it will be 
only necessary to compute the time evolution of one of them, since the others can be found 
by simple permutations. 

The ground state of our model hamiltonian is computed using imaginary-time evolution 
with an initial trial function for <pij(zi, z 2 , z 3 , t = 0) with the required symmetry. The soft 
core potentials parameters were used to fit the energy of the ground state to the experimental 
value E = -7.33 a.u. (199.44 eV), i.e. a = b = 0.4969 (0.262 A) Once the ground state with 
the required accuracy is found, we propagate it in time according to the minimal coupling 
Hamiltonian: 

id 

-QT^aap{zi,Z 2 ,Z 3 ,t) = [H + (pi + p 2 + p 3 )A(t)/c] ® aa p{zu Z 2 , Z 3 , t) (3) 

The vector potential A(t) is assumed to be linearly polarized along the dimension described 
in the model (as usual, the A 2 (t) term of the interaction hamiltonian has been factorized 
as a global phase). As is standard in the Helium case Q], the ionization yield is computed 
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using a partition of the Hilbert space. The extension to the Lithium case reads as: 

Li if \zi | < 15 a.u. and \z2\ < 15 a.u. and | ^3 1 < 15 a.u. 

Li + if \zi\ > 15 a.u. and \zA < 15 a.u. and \zk\ < 15 a.u. , . 

' «. (4) 
Li + if \Zi\ > 15 a.u. and |^-| > 15 a.u. and \zk\ < 15 a .u. 

. Li 3+ elsewhere 

with i, j, fc = 1, 2, 3 

The total ionization yield is obtained by adding the contributions of each of the three 
terms in eq. (0), which describe orthogonal spin configurations. Inspired by the synchrotron 
experiments, we present calculations of the ionization of one-dimensional Li with a pulse 
of frequency of u = 3.016 a.u. (82.06 eV) and intensities ranging from / = 10~ 3 a.u. 
(3.5 x 10 13 W/cm 2 ) up to I = 10 a.u. (3.5 x 10 17 W/cm 2 ). High power coherent radiation 
in this wavelengths are expected to be available at the end of this year in phase 2 of the 
FEL- TTF at Hasylab (Hamburg). To achieve the relevant intensities used in this paper, 
a slight focussing would be needed to focal spots of the order of 10/im. The length of the 
pulse duration (four cycles) is limited by our computer's capabilities. Larger pulse durations 
are expected to increase the ionization yield, but not to alter the fundamental mechanism. 
Note that, specially in the case of shorter pulses, the computations have to be carried out 
over a time interval large enough to allow the ionized population to drift into the proper 
spatial regions (J3|. This interval is typically larger than the pulse length and is determined 
according to the saturation of the ionization yields (see for instance Fig. |2k). 

Figure Q shows the ratio of double to single ionization yields for different intensities, 
computed at about eight laser periods after the end of the interaction. As well established in 
the photoionization of Helium, the change in the slope of this ratio as the intensity increases 
(often referred as knee) is the signature of the transition from non-sequential to sequential 
double ionization. Hence, this figure demonstrates, that double ionization of Lithium also 
shifts from non-sequential to sequential as the intensity increases. In our particular case, 
we may take I t h = 10 15 W/cm 2 as the threshold value between these two mechanisms. 
However, in contrast with the Helium case, in Lithium there are two possible channels of 
correlated double ionization. They correspond to the two different spin configurations of the 
ionized pair: parallel (aa) or antiparallel (af3). Note that the wavefunctions (f>ij(zi, Z2, £3) 
correspond to a definite spin orientation in every coordinate. Therefore, a further partition of 
the spatial volume corresponding to double ionization permits us to track this two channels 
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separately. For instance, in the particular case of (/)i 2 (zi, z 2 , Z3) , the volumes \z±\ < 15 a.u., 
\z 2 \ > 15 a.u., 1 23 1 > 15 a.u., and \zi\ > 15 a.u., \z 2 \ < 15 a.u., \z 3 \ > 15 a.u. describe 
double ionization of an entangled pair with antiparallel spins, while \z%\ > 15 a.u., \z 2 \ > 15 
a.u., 1 < 15 a.u. describes the parallel configuration. Figure |2K shows the dynamics 
of double ionization in each of these two channels, as a function of time at different laser 
intensities. As noted before, the ionized population takes some time to access the spatial 
regions where it is computed. This can be seen in the figure as the ionization yield stabilizes 
at times larger than the interaction. It becomes also apparent the different dynamics of 
ionization for each spin configuration at intensities below the threshold Ith, i.e. when non- 
sequential ionization is the relevant mechanism. On the contrary, both channels tend to 
be equally possible when the ionization is sequential. In conclusion, ionization of electron 
pairs with antiparallel configuration is shown to be more probable when non-sequential 
ionization takes place. Figure |2b plots the relative difference of the ionization yields at the 
end of the computation (final points in Fig. |2K), which is tipically above 50% in the non 
sequential case. This result is in clear agreement with the indication in Q|, in the sense 
that comparison of their experimental results with the ionization of Helium would suggest 
such asymmetry. On the other hand, our previous work in ionization of Helium has shown 
that the e~ — *> 2e~ scattering process is greatly inhibited for the orthohelium case, since the 
parallel spin configuration implies the antisymmetric character of the orbital wavefunction, 
in which Pauli's principle reduces the strength of electron-electron interaction jl^j . Figure 
El demonstrates that this also the case in the double ionization in Li. It shows the density 
distribution corresponding to the term a(l)a(2)/3(3)4>i 2 (zi, z 2 , z 3 ,t) in (J2J) at the end of the 
computation at the planes z\ = 0, z 2 = and z 3 = 0. To improve legibility, black lines outline 
the limits between the regions corresponding to the neutral Li, Li + and Li 2+ . As discussed 
previously, the double ionization is represented by the out-of-axis regions. In this particular 
case, the vertical planes correspond to double ionization of a electron pair with antiparallel 
spins, while the horizontal corresponds to the parallel configuration. The inhibition of this 
later case is, therefore, apparent from this plot. Therefore, and in agreement with jx^J, the 
dominant mechanism of non-sequential double ionization of Lithium at these frequencies 
consists in a first release of an electron, followed almost instantaneously by the scattering 
with one of the two remaining electrons. The exclusion principle makes this scattering most 
effective for the antiparallel spin configuration, hence resulting in a larger ionization yield. 
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In conclusion, we have presented ab initio results for the interaction of Lithium with a 
strong laser field, in a reduced geometry. The model, that has been proven to give deep 
qualitative insight on this process for the simplest atoms (Hydrogen and Helium), is devel- 
oped taking into account the three interacting electrons on equal footing with the proper 
symmetrization of the wavefunction, and including full account of quantum correlations. 
Our results reveal the asymmetry of the non-sequential double ionization process in relation 
with the spin configuration of the entangled ionized pair. We give fundamental insights of 
this phenomena, based in the sensitivity of the electron rescattering to the symmetry of the 
orbital wavefunction. 
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FIG. 1: Double to single ionization ratio. By summing over all spin configurations we can obtain 
the total ionization yield as a function of the intensity. The "knee" structure below the / = 10 15 
W/cm 2 , has been recognized as an indicator of the correlated nature of the ionization process. 

FIG. 2: (a) Double ionization on each of the two possible channels as a function of time. Solid line 
represents the channel with parallel spins and dashed line represents the channel with antiparallel 
spins. For low intensities the yields are different because of the inhibition mechanism, (b) Relative 
difference for each of the channels as a function of the intensity. Channel 1 corresponds to ionization 
with parallel spin | ||) electrons ionized in the final state, and channel 2 correspond to antiparallel 
spin | ||) electrons ionized in the final state. 



FIG. 3: Slices of the three dimensional density of the function \4>n(zi, Z2, z^,t)\ 2 in logarithmic 
scale, corresponding to the spin configuration a(l)a(2)/3(3), at the end of the pulse for I = 10 13 
W/cm 2 . Population in the vertical planes outside the fringes at 15 a.u. (7.9 A ) corresponds to 
double ionization with antiparallel spin electrons (| ||) ). The population in the horizontal plane 
outside the fringes at 15 a.u. corresponds to double ionization with parallel spin electrons (| ||) ). 
Pauli principle inhibits double ionization at the z\ = zi plane, because of the antisymmetry of the 
wavefunction. 
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